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Rabi oscillations under ultrafast excitation of graphene 
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We study coherent nonlinear dynamics of carriers under ultrafast interband excitation of an 
intrinsic graphene. The Rabi oscillations of response appear with increasing of pumping intensity. 
The photoexcited distribution is calculated versus time and energy taking into account the effects 
of energy relaxation and dephasing. Spectral and temporal dependencies of the response on a probe 
radiation (transmission and reflection coefficients) are considered for different pumping intensities 
and the Rabi oscillations versus time and intensity are analyzed. 
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The Rabi oscillations of coherent response under ultrafast excitation of two-level atomic systems were studied 

during last decades, [ll Similar phenomena in bulk semiconductors and quantum wells have also been observed, see 

[3, |3| and reviews. [J, |5| Recently, the properties of graphene under ultrafast interband excitation were investigated, 

and most attention has been concentrated on relaxation dynamics for the case of the linear excitation regime of 

|__i ', epitaxial and exfoliated graphene, see [6|] and |7J respectively. It was found that the relaxation times of photoexcited 

'^ ' electron- hole pairs due to the optical phonon emission are about 0.1 ps [a, 01' so that an investigation of coherent 

r^ ' dynamics is possible during the femtosecond time scales when the regime of Rabi oscillations can be realized with 

' ■ increasing of pumping intensity. To the best of our knowledge, neither experimental nor theoretical treatment of the 

flj I nonlinear coherent response for the Rabi oscillations regime in graphene is not performed yet. Due to the gapless and 

Cj . massless energy spectrum with a neutrinolike dispersion law zLvwP (vw ~ 10^ cm/s is the characteristic velocity for 

. ' the Weyl- Wallace model Q), such a case should be essentially different from the above mentioned cases. ,2-5] 

jrt In this paper, we consider the nonlinear coherent process of carrier photoexcitation in an intrinsic graphene under 

Cj ■ the slow-envelope condition, flTp 3> 1 where J7 is the frequency of light and Tp is the duration of excitation. Eliminating 

I I the nondiagonal components of density matrix, which are responsible for the high-frequency [c>c exp(— if2i)] oscillations 

'■^ ■ of polarization, we describe such a process by the same distribution functions for electrons and holes, fpt because of 

Ch ' the electron-hole symmetry. The kinetic equation for fpt takes form 

o : 

■^' ^^Gifp\t) + Mft\p), (1) 

where G{fp\t) is the interband generation rate and Jo{ft\p) is the collision integral. The general consideration of Eq. 
(1) can be found in 9] and the evaluation of the generation rate for graphene is performed in Ref. 10. Further, we 
solve this equation and analyze the spectral and temporal dependencies of the transmission and reflection coefficients 

(y-^ for a probe radiation at different pumping levels. 

fS| ' The coherent interband photoexcitation caused by the in-plane electric field WtE exp(— ifit) -I- c.c. with the envelope 

. -/ I form- factor wt is described by the generation rate in Eq. (1) [l^] 



X cos 



f^-n]t' 



(1-2/pt+r). (2) 



Here we neglect the multi-photon interband transitions under the condition {eEvw /fiQ?)"^ ^ 1 ^ind the dephasing 
relaxation time Td is introduced phenomenologically. The non-local factor (1 — 2fpt+t') describes the Pauli blocking 
and the temporal memory effect which is responsible for the Rabi oscillations. Since the cascade emission of optical 
phonons dominates in the energy relaxation, [6, 7, 11] Eq. (1) involves the collision integral 

■h{It\p) — ^'p+po(l ~ fpt)fp+pot — ^'p_p(,(l — fp-pat)fpt, (3) 

where po = Tiujq/vw is the characteristic momentum, Tiloq is the optical phonon energy, and Vp — VQp/fi is the 
relaxation rate for the spontaneous emission of optical phonons. Under the condition Vp^Tp <^ 1 we neglect interband 
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transitions due to optical phonon emission. We also restrict ourselves by the model with a single phonon of energy 
TiojQ ~0.2 eV with the efficiency of coupling determined by the characteristic velocity Vq ~ 10^ cm/s. These parameters 
are in agreement with the results of calculations of the electron-phonon interaction [1J| and with the measurements of 
relaxation dynamics. [aQ The problem given by Eqs. (1) - (3) is solved below with the initial condition /pt-y-oo = 
which is correspondent to the undoped graphene. 

The analytical solution of the problem formulated can be found for the coUisionless case, r^ — > oo and vq -^ 0, 
under the resonant condition p — >■ hfl/2v\Y = Po, when the oscillating factor is absent in the integral generation rate 
(2). For such a case, the integro-differential Eq. (1) can be transformed into the second order differential equation 
and the solution takes form 



U- 



'■pn,t 



cos 



f77l— I dt' 

V2/ex / Wt' 



(4) 



Here we introduced the dimensionless intensity, /e; 



(eETpVw/^^)'^, so that the Rabi oscillations of the resonant 
E) takes place. At i — >■ oo one obtains the resonant 



distribution with time and with field strength (y/2Tf^x oc 

distribution /p=pn,t-i-oo = (1 ^ cos.4ea;) /2 which is determined by the dimensionless area of the incident pulse A 

V2IexJ^^dtWt/Tp. 
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FIG. 1: Contour plots of photoexcited distribution fpt versus dimensionless momentum and time, (p — pii)/Ap and t/rp for 
the coUisionless regime at pumping levels corresponding to Aex = ti" (a), 271" (b) and Stt (c). 



The coUisionless case at p y^ pn is described by the reduced equation (1) dfpt/dt — G{fp\t). The numerical solution 
of this equation is obtained here with the use of the finite difference method [1J| and the Gaussian form-factor 



141 In Fig. 1 we plot fpt versus dimensionless time, t/Tp, and momentum {p-'Pfi)/Ap 
Here Ap — h/(vwTp) determines a width of photoexcited distribution at i > 2Tp while a 



wt = \/2/tt exp - (t/TpY 

which is centered at p^. 

width of distribution at i ~ increases with pumping intensity, as it is shown in Figs. 1 a-c. As lex increases, a 

temporal Rabi oscillations dX p ^ pn are similar to the transient evolution described by Eq. (4). 

We turn now to consideration of the problem (1) - (3) taking into account the dephasing and energy relaxation 
processes, when r^ and uq are finite. In the calculations below we use the pumping frequency TiU, =1.5 eV and the 



dephasing time r^ 



=85 fs, which is correspondent to the process of spontaneous emission of optical phonons. 



Since the cascade emission of dispersionless optical phonons, described by Eq. (3), a multipeak structure of fpt with 
maxima centered around pn — kp^ (fc = 0, 1, . . .) is realized. The only first phonon repetition is essential during the 
photoexcitation process, when t < 2Tp, see Figs. 2a and 2b where the contour plots of the photoexcited distributions 
fpt at the pumping level Aex = 37r are shown for the pulse durations Tp =30 fs and 80 fs. One can see that the 
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FIG. 2; Contour plots of photoexcited distributions fpt versus energy vwp and time for pulse durations Tp =30 fs (a) and 80 
fs (b) at pumping level corresponding to Aex = Stt. 

photoexcited distribution appears to be narrower with increasing of Tp and that the Rabi oscillations are visible at 
Td ~ Tp (Fig. 2b). 
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FIG. 3: (Color online) Transient evolution of concentration of photoexcited carriers for pulse durations Tp =30 fs (a) and 80 fs 
(b) at pumping levels corresponding to Aex = tt (1), 2% (2), 37r (3), An (4), Stt (5), and Bvr (6). (c) Photoexcited concentration 
versus Aex at t ^ Tp for different Tp (marked). 



Using the distribution fpt obtained, we calculate here the concentration of photoexcited carriers 



nt = -^72 / dppfpt 

TTfi Jo 



(5) 



versus time and versus pumping level at t ^ Tp. As one can see from Figs. 3a and 3b, the amplitudes of temporal 
oscillations decreases with increasing of Tp but the oscillations remain visible at Tp ^ tj_. The oscillatory behavior of 
nt^oo versus pumping level (which is proportional to A^^) appears to be suppressed at Tp ^ tj, as it is shown in Fig. 
3c. The condition Aex = tt corresponds to the pulse energies ^15 nJ or ~5.6 nJ for Tp =30 fs or 80 fs and for the 



sport area ~ 10 cm (note, that the pulse energy oc t and 



if A 



mr). 
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FIG. 4: (Color online) Spectral dependencies of the dynamic conductivity, Keauit (a) and ImAauit (b), around the pumping 
energy 1.5 eV for delay times between —30 and 40 fs with step 10 fs (marked by 1 - 8) at pumping level corresponding to 
Aex = Stt and pulse duration Tp =30 fs. 

The transient response on a probe radiation ex exp(— zwi) is determined by the dynamic conductivity cr^t which 
was evaluated in ^9, 15] for the collisionless case ojTd 3> 1, when the parametric dependency on time takes place. The 
real and imaginary parts of a^t are given by 



Rea., = -(1 



Imcr^t = (Tuj ^V 



nh 






(6) 
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where V means the principal value of integral and we introduced the time-independent contribution a^, described the 
undoped graphene in the absence of photoexcitation (below we use the phenomenological expression for W^ introduced 
in [ig). The negative absorption condition Kea^t < takes place if fp^.t > 1/2; these regions were marked by thick 
curves in Figs. 1 and 2. The peak shape of Kea^jt centered at w = fJ leads to a visible spectral dispersion of Ima^t 
at w < 17 and lo > fl. 

The reflection and transmission coefficients, Rc^t and T^^t, are written through the dynamic conductivity as follows: 

M 



Re 



1 - Ve- 47rcr^t/c 
1 + \/f^ + 4:nauit/c 

1 + Ve + 47rCT^t/c 



(7) 



Tujt — 



Here e is the dielectric permittivity of a thick substrate and the geometry of normal propagation of radiation was 
considered. The contour plots of R^jt and T^jt versus t and lo are shown in Fig. 5 at pumping level corresponding to 
Aex = 3tt for Tp = 30 fs and 80 fs. Temporal oscillations of response correlate with the evolution of a^t (Fig- 3) due 
to oscillations of distribution shown in Fig. 2. Spectral width of Rabi oscillations appears to be broader for Tp — 30 
fs (c.f. Figs. 5a, 5b and 5c, 5d with Figs. 2a and 2b). The oscillatory behaviors of the reflection and transmission 
coefficients at the pumping frequency (oj = i7), R^t and Tnt, versus pumping level for the delay times Tp, 2rp, and 
3Tp are shown in Figs. 6a and 6b, respectively. One can see a few-percent modulation of response at t > Tp versus 
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FIG. 5: Contour plots of reflection (a, c) and transmission (b, d) coefficients versus energy of probe radiation huj and time for 
pulse durations Tp =30 fs (a, b) and 80 fs (c, d) at pumping level corresponding to Aex ~ Stt. 
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FIG. 6: (Color online) Reflection (a) and transmission (b) versus intensity for Tp = 30 fs and 80 fs (solid and dotted curves) at 
different delay times t = Tp, 2rp, and 3rp (marked) and for u = f2. 



pumping intensity. A visible damping of this modulation takes place due to the cascade emission of optical phonons, 

Our calculations are based on the follovifing assumptions. Since the Coulomb renormalization of interband transitions 
is weak, [l7[ we have used a single-particle approach. All homogeneous dephasing mechanisms, including optical 
phonon emission and carrier-carrier scattering, have been described phenomenologically, through the dephasing time 
Td in Eq.(2). An inhomogeneous broadening due to long-scale disorder is not taken into account, so that the results 
are valid for a high-quality graphene. The model collision integral (3), which is written through an effective phonon 
energy tiojo and a relaxation frequency Up, is used because the only first step of cascade emission of phonons is essential 
during the Rabi oscillation process. A detailed calculations can improve a precision of vq used in Eq. (3) and will give 
a widened peak of distribution at pn — Pq- The above-listed simplifications of the damping processes do not change 
the temporal dynamics under consideration and we have demonstrated that the Rabi oscillations are observable in 
a typical graphene sample. The rest of assumptions [the geometry of normal propagation of radiation, the Gaussian 
shape of excitation, and the collisionless approximation used in Eqs. (6)] are rather standard. 

Summarizing, we have described the mechanisms of coherent nonlinear response of an intrinsic graphene under 
ultrafast interband excitation. The results obtained demonstrate that the Rabi oscillations, both versus time and 
versus pumping intensity, can be easily observed for femtosecond time scales (up to 0.1 - 0.2 ps) at pumping intensities 



■^3 - 30 GW/cm^ (which are correspondent to pulse energies '-^lO - 100 nJ for the sport area ~ 10 
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